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\Q . Abstract 
\& • A model of holographic dark energy with an interaction with matter fields has been investigated. Choosing 

o 

' the future event horizon as an IR cutoff, we have shown that the ratio of energy densities can vary with 

time. With the interaction between the two different constituents of the universe, we observed the evolution 
of the universe, from early deceleration to late time acceleration. In addition, we have found that such an 
interacting dark energy model can accommodate a transition of the dark energy from a normal state where 
wd > —1 to wd < —1 phantom regimes. Implications of interacting dark energy model for the observation 
of dark energy transition has been discussed. 
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The present acceleration of the universe expansion has been well established through numerous 
and complementary cosmological observations A consistent picture has indicated that nearly 
three quarters of our universe consists of a "dark energy" , which is responsible for the accelerated 
expansion. However, the nature of such a dark energy is still rather uncertain. Explanations have 
been sought within a wide range of physical phenomena, including a cosmological constant, exotic 
fields, a new form of the gravitational equation, new geometric structures of spacetime, etc Q, Q|< 
Recently, a new model stimulated by the holographic principle has been put forward to explain 
the dark energy J] . According to the holographic principle, the number of degrees of freedom of a 
physical system scales with the area 

in quantum field theory a short distance cutoff is related to a long distance cutoff due to the limit 
set by formation of a black hole, which results in an upper bound on the zero-point energy density. 
In line with this suggestion, Hsu and Li jfjjQ] argued that this energy density could be viewed as 
the holographic dark energy density satisfying po = 3c 2 /L 2 , where c 2 is a constant, L is an IR 
cutoff in units M 2 = 1. Li [J] discussed three choices for the length scale L which is supposed to 
provide the IR cutoff, such as the Hubble radius, the particle horizon and the event horizon. He 
demonstrated that only identifying L with the radius of the future event horizon, we can get the 
correct equation of state of dark energy and obtain the desired accelerating universe. Besides other 
applications of the holographic principle in cosmology 7], the holographic dark energy model is a 
new example showing that holography is an effective way to investigate cosmology. Related works 
on the holographic dark energy can be found in 8]. The holographic dark energy model is found 
in consistent with the observational data 

Available models of dark energy differ in the equation of state parameter wd as well as the 
variation of wd itself during the evolution of the universe. The cosmological constant with wd = 
— 1, is located at a central position among dark energy models both in theoretical investigation 
and in data analysis jLO]. In quintessence, K-essence, Chaplygin gas and holographic dark energy 
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models |2j, |4J, wd stays 



energy have wd < 



jigger than —1, and cannot cross —1. The phantom models of dark 



111 ]. Recently, the analysis of the type la supernova data indicates that 
the time varying dark energy gives a better fit than a cosmological constant 12]. These analysis 
mildly favor the evolution of the dark energy parameter wd from wd > —1 to < —1 at recent 
stag e. Although the galaxy cluster gas mass fraction data do not support the time- varying wd 
13], theoretical attempts toward the understanding of the wd corssing —1 phenomenon have been 
started. Some dark e nerg y models, such as the one containing a negative kinetic scalar field and 



a normal scalar field 



141 ] . or a single scalar field model |l5j], a Gauss-Bonnet brane world with 



2 



induced gravity 16] and the generalized implicitly denned dark energy equation of state model 



[13| , have been constructed to gain insight into the occurrence of the transition of the dark energy 
equation of state and the mechanism behind this transition. Other studies on the wd crossing —1 
can be found in jisl |. 

Most discussions on dark energy rely on the fact that its evolution is independent of other 
matter fields. Given the unknown nature of both dark energy (DE) and dark matter (DM), which 
are two major contents of the universe, one might argue that an entirely independent behavior of 
DE is very special [ljj- 
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special [19|]. Studies about the interaction between DE and DM have been carried out 



21, Q. It was argued that the interaction will influence the perturbation dynamics 



and could be observable through the lowest multi-poles of the CMB spectrum [20|]. Recently, by 
considering the interaction between DE and DM in the holographic DE model, it was even found 
that the Hubble scale might be used as the IR cutoff to explain the acceleration of our universe 
j^ . In this work, we are going to extend the inclusion of interaction between DE and DM into 
the holographic DE model with the future event horizon as the IR cutoff. The ratio of energy 
densities can be varied with time. As a result, we find that with the interaction between DE and 
DM, the model can give an early deceleration and a late time acceleration. In addition, in this 
holographic model, the appropriate coupling between DE and DM accommodates the transition of 
the DE equation of state from wn > — 1 to wd < — 1- This property could serve as an observable 
feature of the interaction between DE and DM in addition to its influence on the small I CMB 
spectrum argued in j^ . 

The total energy density is p = p m + pp, where p m is the energy density of matter and pn = 
3c 2 R E 2 is the dark energy density. Here, the Planck mass is taken as unit. We followed |J{ by 
choosing the future event horizon, Re = a j^z, as the IR cutoff, where c 2 is a constant. The 
total energy density satisfies a conservation law. However since we consider the interaction between 
DE and DM, p m and po do not satisfy independent conservation laws, they instead satisfy 

p m + 3Hp m = Q (1) 

and 

p D + 3H(1 + w D )p D = -Q , (2) 

where wd is the equation of state of DE, Q denotes the interaction term and can be taken as 
Q = 3b 2 Hp with b 2 the coupling constant 22]. This expression for the interaction term was 
first introduced in the study of the suitable coupling between a quintessence scalar field and a 



pressureless cold dark matter field [l^]. The choice of the interaction between both components 
was to get a scaling solution to the coincidence problem such that the universe approaches a 
stationary stage in which the ratio of dark energy and dark matter becomes a constant. In the 
context of holographic DE model, this form of interaction was derived from the choice of Hubble 
scale as the IR cutoff 22 1. 



Taking the ratio of energy densities as r = p m /pD, from © and © we have 

r = 3b 2 H(l + r) 2 + 3Hrw D . (3) 



In |22|, the IR cutoff was chosen as being the Hubble scale, which leads to a constant r. We 
choose the future event horizon as the IR cutoff so that r is no longer a constant in the holographic 
DE model. 

Using the Friedmann equation, O m + f2£) = 1, where Q m = p m /(3H 2 ) and = pd/(3H 2 ), we 
have r = (1 — Qd)/Q,d and r = —ClD/Qfj- Combining with (j3J), we get the equation of state of 
DE, 

w D = [-h D /n 2 D -3b 2 H(l+r) 2 }/(3Hr) 

= I (4) 

where the dot is the derivative with respect to time and the prime is the derivative with respect 
to x = In a. 



From the Friedmann equation, the future event horizon can be expressed as Re = c\/l + r/H = 
a J* t °° dt/a. Taking the derivative with respect to t and using the definitions of r and r as well as 
equation ((3}, we arrive at 



tt' D „ „ Nr 1 2 3b' 
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Back to (JIJ) we get now 

WD = -1/3 - 2y / fW(3c) - b 2 /Q D . (6) 

Neglecting the interaction between DE and DM, © is equivalent to equation (18) of Li's paper in 
J]. We will show that the interaction between DE and DM brings rich physics. 
The deceleration parameter can be expressed in this model as 

q = -da/a 2 = -H/H 2 - 1 

= 1/2-36 2 /2-^d/2-^ /2 /c , (7) 



where H/H 2 = 3(-l -w D - r)Q D /2. 

With eqs (jlEJ), we are in a position to investigate the evolution of the DE and its influence on 
the expansion of the universe. 

Since the DE plays a more important role in the evolution of the universe with the flow of 
cosmological time, we require £l' D > 0, which leads to 
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Figure 1: Evolution of the DE for a fixed interaction parameter with DM (b 2 ) but for different values of the 
constant c. 

The behavior of the DE evolution puts an upper limit on the interaction between DE and DM, 
which has not been observed in the study of the holographic DE model by taking the Hubble 
scale as IR cutoff decreases with the increase of c. Choosing two different values of c 

with ci < C2, their corresponding largest allowed values b 2 ^^ are b\_ max and b\_ max and satisfy 
bi-max > ^2-max- Taking b 2 as a common allowed value of the coupling for c\ and C2 cases, it is 
easy to see that b 2 s is closer to b\_ max than b 2 _ max , which means that in the large c case the same 
allowed coupling b 2 is stronger than in the small c case. 

The dependences of the evolution of DE with respect to the constants b 2 , c are shown in Fig. f 
and Fig. 2, respectively. From Fig. 1 we learn that for the fixed coupling between DE and DM 
within the allowed range (JHJ), the DE starts to be effective earlier when c is larger. Since the same 
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Figure 2: Evolution of the DE with a fixed constant c but different values for the coupling with DM. 

b 2 corresponds to a stronger coupling between DE and DM for bigger c cases, £Id will tend to a 
smaller value at late stage when c is bigger due to this stronger coupling. 

For a fixed c, the dependence of on b 2 is shown in Fig. 2. When b 2 is larger, at early stage, 
the DE starts to be effective earlier. However at a later stage, since the coupling between DE and 
DM becomes stronger, Ud approaches a smaller value for a larger b 2 . 

Including the interaction, our model naturally shows that our universe has an accelerated ex- 
pansion in the late stage and on the other hand it also displays a deceleration in the early era. In 
Fig. 3 we show that for the same coupling between DE and DM, the acceleration starts earlier for 
larger c, since the DE develops earlier for larger c with the same coupling b 2 . In Fig. 4 we show 
that for the same c, the acceleration starts earlier for larger b 2 . This is also due to the fact that 
for the same c, DE develops earlier for bigger b 2 as shown in Fig. 2. 

We now discuss the equation of state of the DE with the interaction between DE and DM. From 
we learnt that wd has a maximum value -1/3 - (b/V3c) 2 / 3 at n D - cr = (36 2 c) 2 / 3 . 

To allow the DE transition as indicated by recent observations with wd crossing the border —1 and 
wd < —1 at present stage, we need b 2 > b 2 r = 2Qd(1 — \/TTd/c)/3, where b 2 r is got from wd = — 1- 
Such b 2 r should be in the reasonable range (jBJ). Thus we get c < 2V^d/(3^_d — 1). Meanwhile b 2 r 
should be positive, which leads to c > \ATd- According to the holographic principle, the entropy 
of the universe is bounded by S = ttR e , where Re = c/ (H^/Qd)- If we require the entropy of the 
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Figure 3: Dependence of the deceleration parameter on the constant c for a fixed coupling between DE and 
DM. 




Figure 4: Dependence of the deceleration parameter on the coupling between DE and DM for a fixed 
constant c. 

universe do not decrease, we need Re = c/V^d — 1 > 0. thus we get c > V^D- Therefore the 
lower bound on c is exactly the requirement of the second law of thermodynamics discussed in jj]. 
Since the DE evolves independently of DM in cases (b 2 = 0), wd < — 1 requires c < \/Qd, which 
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violates the second law. Thus in the holographic model with entirely indepedent behavior between 
DE and DM, it is impossible to have wd crossing —1 allowing the phantom energy to exist in the 
late stage of the universe. 

Further examining the function b 2 r , we see that it has a maximum value b 2 r _ max = 8c 2 /81 at 
f&D = 4c 2 /9. In order to have wd < —1, we need b 2 > b 2 r . However the coupling b 2 cannot be 
arbitarily large, since when b 2 > b 2 r _ max and wd < — 1 always. The upper bound on b 2 can also 
be got by requiring wo-max > — 1 thus allowing wd to cross —1 and to stay below —1 later. In 
order to have wjj < —1 with Qdo = 0.7, such a maximum value b 2 r _ max should appear early when 

= 4c 2 /9 < 0.7. Combining constraints on c, we have 



fl D <c< 1.255 . (9) 

Now let us focus on the parameter space of b 2 , which we have already obtained, b 2 r < b 2 < 
8c 2 /81. Meanwhile taking account of the behavior of the wn function, we also require f2_D_ cr < ^_do> 
which leads to b 2 < 0^q/(3c), allowing the wd transition from wd > —1 to wd < —1 in recent 
times. Using Qdo = 0.7, the combined constraint on b 2 reads 

1.4(1 - VoJ/c)/3 < b 2 < 8c 2 /81 , (10) 

to accommodate wd crossing —1. 

For a fixed c, the dependence of wd on b 2 is shown in Fig. 5. We see that for larger b 2 , wd 
crosses —1 earlier. For a too small b 2 , wd crosses —1 too late. Actually these small b 2 are over the 
lower bound of (10) and are not consistent with the observation. 

From the future precise observation of the location of the transition of the wu from wd > —1 
to wd < — 1, it is possible to understand the interaction between DE and DM. In the future this 
could serve as another observable feature of the interaction between DE and DM in addition to 
the low CMB spectrum discussed in j^ . 

We have also fitted our model with the golden SN data. We got Q m0 = 0.39+°A6' 1)2 = 0- 00 -aoO' 
c = OAOtoll and x 2 = 174.44. If we fix c = 1, we have O m0 = 0.26±g;^|, b 2 = O.OOl^, and 
X 2 = 177.08. This shows that our model is consistent with the SN data and the constraint on 
parameter spaces we discussed are compatible with the observations. 

In summary, we have studied the holographic dark energy model with an interaction between 
DE and DM. Choosing the future event horizon as the IR cutoff, we obtained a generalized model 
with the time dependent ratio of energy densities which cannot be realized just by adopting the 
Hubble scale as the IR cutoff. With the interaction between DE and DM, a richer physics arises. 
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Figure 5: Behavior of the equation of state of DE for different couplings between DE and DM. 

In addition to showing the comprehensive history of the evolution of our universe from the early 
deceleration to late acceleration, we have also found that the interaction between DE and DM 
can accommodate the transition of the equation of state of DE from wd > — 1 to i«d < — 1. We 

have constrained the parameter space of our model to explain the observations. We argued that 

l I 

in addition to some observational feature for small I CMB spectrum |20j, the interaction between 
DE and DM could also be observed by the wd crossing —1 behavior in the future. Comparison 
to the golden SN data has been made, and we showed that our model is also consistent with such 
observations. 
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